Fermi-surface evolution across the magnetic phase transition in the Kondo lattice 

model 
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We derive, by means of an extended Gutzwiller wavefunction and witliin the Gutzwiller approx- 
imation, the phase diagram of the Kondo lattice model. We find that generically, namely in the 
absence of nesting, the model displays an /-electron Mott localization accompanied by a discon- 
tinuous change of the conduction electron Fermi surface as well as by magnetism. When the non 
interacting Fermi surface is close to nesting, the Mott localization disentangles from the onset of 
magnetism. First the paramagnetic heavy fermion metal turns continuously into an itinerant mag- 
net - the Fermi surface evolves smoothly across the transition - and afterwards Mott localization 
intervenes with a discontinuous rearrangement of the Fermi surface. We find that the /-electron 
localization remains even if magnetism is prevented, and is still accompanied by a sharp transfer 
of spectral weigth at the Fermi energy within the Brillouin zone. We further show that the Mott 
localization can be also induced by an external magnetic field, in which case it occurs concomitantly 
with a metamagnetic transition. 

PACS numbers: 74.20.Mn, 71.27.-|-a, 71.30.+h, 71.10.Hf 



I. INTRODUCTION 

The physics of heavy-fermion compounds remains a 
fascinating and challenging issue within strongly corre- 
lated materials. Recently, considerable experimental and 
theoretical interest has focused on the physical behav- 
ior across the magnetic quantum phase transition that 
is traditionally expected to occur when the Ruderman- 
Kittel-Kasuya-Yosida (RKKY) interaction overwhelms 
Kondo screening.^ This transition is induced experimen- 
tally by external parameters like chemical composition, 
pressure or magnetic field, see for instance Refs. d andU 
as well as references therein, and is commonly accompa- 
nied by topological changes of the Fermi surface^'^'^'^ 
and anomalous behavior of various transport and ther- 
modynamic quantities. ^•■^ The theoretical debate on this 
subject has so far mainly followed two different direc- 
tions.-^' One ascribes the changes of the Fermi surface to 
an /-electron Mott localization,^ which is assumed to oc- 
cur concomitantly with magnetism as well as with the 
appearance of transport and thermodynamics anoma- 
lieSf^i^'^i^ The alternative proposal assumes that mag- 
netism is predominantly an instability of an itinerant 
phassi^ii^, hence that the Fermi surface changes arise 
simply by the spin polarization of dispersing band o^^i^^ 
and the anomalous behavior by critical magnetic quan- 
tum fluctuationSflS This issue has been very recently ad- 
dressed theoretically in the periodic Anderson model by 
De Leo, Civelli and Kotliar-iii^ using a cluster exten- 
sion of dynamical mean field theory (CDMFT). Upon 
decreasing the hybridization between /-orbitals and con- 
duction electrons, a weak first order phase transition from 
a heavy-fermion paramagnet to an itinerant antifcrro- 
magnet has been found. Remarkably, when these au- 
thors force CDMFT not to break spin SU(2) symmetry 



and follow the metastable paramagnetic solution, they 
find an orbital-selective Mott localization - a pseudogap 
opens in the /-electron spectral function at the chemical 
potential, although low energy spectral weight remains 
within the Mott-Hubbard gapi^ii^ - for a hybridization 
between / and conduction electrons almost coincident 
with the value at which, allowing for magnetism, the an- 
tiferromagnetic transition occurs. This result suggests 
that the magnetic phase transition masks an incipient 
Mott localization of the /-electrons, which could be- 
come visible above the Neel temperature or by suppress- 
ing antiferromagnctism. A complementary attempt has 
been almost contemporarily performed by Watanabe and 
Ogata^i These authors analyse by a variational Monte- 
Carlo (VMC) technique a Gutzwiller wavefunction for a 
Kondo lattice model in a two-dimensional square lattice. 
The variational phase diagram as function of the Kondo 
exchange depends non-trivially on the electron density. 
Very close to the compensated regime (one conduction 
electron per impurity-spin), upon decreasing the Kondo 
exchange there is first a second-order paramagnetic-to- 
antiferromagnetic phase transition, followed by a first- 
order transition between two magnetic phases with dif- 
ferent Fermi surfaces. Moving away from the compen- 
sated regime, the second order phase transition disap- 
pears and they find a single first order line separating a 
paramagnetic phase from an antiferromagnetic one with 
different Fermi surfaces. These VMC results suggest that 
the magnetic transition and the topological change of the 
Fermi surface are not necessarily coincident, which has 
been also observed in a very recent experiment. Since a 
variational calculation can only access ground state prop- 
erties and not subtle dynamical features like an orbital- 
selective Mott transition, and keeping into account the 
differences between the periodic Anderson model and the 



2 



Kondo lattice model, the VMC21 and CDMFTl^; results 
might not be incompatible one to the other, and instead 
describe the same physical scenario although from two 
different perspectives. Should this be the case, it would 
undoubtly represent a step forward in the comprehension 
of heavy-fermion physics. To settle this question, one 
should for instance try to get closer to the compensated 
regime by CDMFT and check whether the /-localization 
and the on-set of magnetism disentangle from each other 
as predicted by VMC. Alternatively, one could carry on 
with variational calculations trying to uncover features 
that indirectly signal the /-localization. This is the aim 
of the present work. We note, by the way, that finite av- 
erage values of the hybridization between / orbitals and 
conduction electrons, in the periodic Anderson model, 
or of the Kondo exchange, in the Kondo lattice model, 
must not be interpreted as absence of /-localization in a 
proper variational calculation, since the hybridization or 
the Kondo exchange are part of the Hamiltonian. There- 
fore other quantities must be identified that are accessible 
by a variational calculations. 



In particular, in this work we adopt a variational tech- 
nique based on a recent multi-band extension^"^ of the 
so-called Gutzwiller approximation to evaluate analyti- 
cally average values on Gutzwiller variational wavefunc- 
tions.^*^'-^ This method is not exact like VMC, unless in 
the case of infinite-coordination lattices. However, we 
have found that a variational wavefunction richer than 
that of Ref. Hf] seems to compensate for the approxima- 
tion adopted to calculate average values, thus leading to 
the same phase-diagram as the one obtained by VMC 
in the case of a two-dimensional square lattice^i. En- 
couraged by this result, we have extended the analysis of 
Ref. [2l|. Specifically, we have derived the phase diagram 
forcing the wavefunction to remain paramagnetic. Simi- 
larly to the CDMFT calculation of Ref.[ll we have found 
that a transition accompanied by a topological change 
of the Fermi surface exists also in this case, although 
is masked by magnetism when we allow for it. Finally, 
we have analyzed the role of a uniform magnetic field in 
the paramagnetic phase and found a matamagnetic in- 
stability near the above phase transition, suggestive of 
the metamagnetism observed experimentallyi^ 



The plan of the paper is as follows. In section II we 
introduce the model and the variational technique. The 
variational phase diagram is presented in section III. In 
section IV we discuss the properties of the variational 
wavefunction in the paramagnetic sector, while in section 
V we consider the effect of a magnetic field in connection 
with metamagnetism. Section VI is devoted to conclud- 
ing remarks. Finally, in the appendix we present some 
technical details of the variational method employed. 



II. THE MODEL AND THE VARIATIONAL 
METHOD 

We consider a Kondo lattice model (KLM) described 
the Hamiltonian 

<RR'> (T ^ ^ 

S/R-S,R = 7^o+Hj, (1) 

R 

where creates a conduction electron at site R with 
spin a that can hop with amplitude —t to nearest neigh- 
bor sites, S/R is the spin-1/2 operator of the /-orbital 
and ScR the conduction electron spin-density at site R. 
In what follows, we assume a bipartite lattice. To study 
this Hamiltonian we introduce the following variational 
Gutzwiller wavefunction 

R 

where |^'o) is the ground state of a non-interacting two- 
band variational Hamiltonian describing hybridized c and 
/ orbitals, while Vr is a local operator that modifies the 
relative weights of the local electronic configurations with 
respect to the uncorrelated wavefunction. In particular, 
we will assume for V-r. the general expression 

7'r = 5I ^r„(R)|r,R)(n,R|, (3) 

r,n 

where |r, R) and |ri, R) span all electronic configurations 
of the c and / orbitals at site R, with the constraint that 
the states |r,R), but not |n, R), have just a single /- 
electron. 

The variational wavefunction ^ has been widely used 
to study the periodic Anderson model as well as its strong 
coupling counterpart, the Kondo lattice model, within 
the Gutzwiller approximationj ^^i^^i^^i'^°i^^i'^^ However, in 
all the earlier works the operator Pr has been choosen 
to act only on the /-orbitals states. For instance, in 
the KLM that we consider, this choice would reduce to 
take Vr as the projector onto singly occupied /-orbitals, 
namely 

7'R-^|^,R)(^,R| = (n/RT-n/Rx)^ (4) 
r 

where ny^Ro- = /ro-Zro-- This assumption implies that 
the spin correlations induced by the exchange J in ([T|) 
are only provided by the uncorrelated wavefunction I'I'o). 
The more general form of Pr, Eq. ([3]), that we assume in 
what follows, permits to include additional correlations 
besides those included in the wavefunction |^o)) in par- 
ticular the tendency of the conduction electrons to couple 
into a singlet with the localized spins. 

The variational procedure amounts to optimize both 
the parameters Arn(R) as well as those that identify \'i'o) 
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by minimizing the average value of the Hamiltonian ([T]). 
In general this task can be accomplished only numeri- 
cally, for instance by means of VMC as actually done 
by Watanabe and Ogata^^ with the simple choice of Vu 
as in Eq. (|4]). However, in infinite coordination lattices 
many simplications intervene that allow to evaluate av- 
erage values analyticallyi^i^i^ In this work we follow 
an extension^'^ of the multi-band method developed by 
Biinemann, Weber and Gebhard^^i^^ that allows to han- 
dle with non-hermitean operators 'Pr,, which is generally 
the case since the bra (n, R| in ([3]) can have any number 
of /-electrons while the ket |r, R) is forced to have only 
one. 

We start assuming that Vr is not the most general as 
possible but is subject to the following two conditions^! 
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l-'Rcr — I pf 



Rct'^Rct /rct/rct 



(5) 
(6) 



(7) 



is the local single-particle density matrix operator. If 
Eqs. ([5]) and ([6]) are satisfied, then one can show^^^i^ 
that, in an infinite-coordination lattice, the average value 
of ^ that has to be minimized is 



E 



<RR'>cr 

KcA^') cn'. + /rv ) + H.c. 

^E(*ol^RS/R.SeR7'j,|*o; 
R 



Rcr 



(8) 



The hopping renormalization coefficients Z are obtained 
through the following equations: 



(*0|^^cL^r/rJ*0> 



CrctI^o) 

Z,,,(R) (^oIcL/rJ^o) + ^c/.(R) (*o|/L/rJ*o> 

I 



(9) 
(10) 



Therefore the variational calculation reduces, in infinite 
coordination lattices and provided Eqs. ^ and ^ are 
satisfied, to calculate expectation values on the Slater 
determinant uncorrelated wavefunction, which is analyt- 
ically feasible since Wick's theorem applies. In the ap- 
pendix we show how one can manipulate and simplify 
all the above expressions to get manageable formulas. 
For finite-coordination lattices the expression ([5]) for the 
variational energy is nomore correct. However, it is com- 
mon to keep using also in these cases the same formula 
([5]) with the coefficients Z defined above - the so-called 
Gutzwiller approximation^liS^. Moreover, it is also com- 
mon to intepre t^^'^^ the non-interacting Hamiltonian , 
see Eq. (|A.2ip in the appendix, whose ground state is the 
optimized wavefunction |^'o), as the Hamiltonian of the 
quasiparticles within a Landau-Fermi liquid framework. 

Before moving to the presentation of our variational re- 
sults, we want to mention some important consequences 
of choosing that acts both on the / and on the 
c orbitals. A drawback of the conventional Gutzwiller 
wavefunction with of Eq. ([5]), which was pointed 
out already by Fazekas and Miiller-Hartmann in Ref. Is^ . 
is that, for small J, the paramagnetic solution gains a 
singlet-condensation energy that has a Kondo-like ex- 
pression oc exp(— 1/Jp), with p the conduction electron 
density of states at the chemical potential. On the con- 



trary, any magnetic solution gains a local exchange en- 
ergy of order J^p - the average value of ^X^r ^/R. ' ^cR 
- independently of the spatial arrangement of the mag- 
netic ordering. This result would remain true even for a 
single impurity Kondo model and is obviously incorrect. 
Our wavefunction partially cures this deficiency because 
Vr is able to induce additional spin-correlations among 
c and / electrons, although only locally. 



We further note from ^ that the action of the 
Gutzwiller operator effectively generates an intersite 
hopping between the /-electrons, absent in the original 
Hamiltonian ([1]), which correlates different sites hence 
can play an important role in determining the topology 
of the Fermi surface as well as in stabilizing magnetic 
structures. Even though our method for computing aver- 
age values is not exact in finite-coordination lattices, the 
more involved form of Vj^ of Eq. ^ with respect to ^ 
partly compensates for this weakness - the variational 
Hamiltonian contains inter-site /-/ and f-c hopping - 
leading to results that are very similar to those obtained 
by exact VMC, as we are going to show. 
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FIG. 1: (Color online) Variational phase diagram as function 
of the conduction electron density ric and of the Kondo ex- 
change in units of half the bandwidth, J/D. The solid line 
with the circles represents a first order line, while the dot- 
ted line is a second order transition. The error bars along 
the second order phase transition line reflect the variational 
uncertainty of a precise location of the continuous transition. 
The same problem does not arise along the discontinuous first 
order line. PM stands for paramagnetic heavy-fermion metal, 
while AF stands for an itinerant antiferromagnet , the sub- 
scripts "e" and "h" are borrowed from Ref. [2l| and refer to 
the electron- like, "e", or hole-like, "h", character of the Fermi 
surface, see Fig. [3] 



III. VARIATIONAL PHASE DIAGRAM 



In Fig. [T] we show the variational phase diagram as 
function of the Kondo exchange J, in units of D, versus 
the conduction electron density < ?ic < 1- Close to the 
compensated regime Uc = I, one conduction electron per 
spin, we do find, similarly to Watanabe and Ogata^^, two 
successive transitions as J/D is reduced from the heavy- 
fermion paramagnetic phase. First, Neel antiferromag- 
netism appears by a second order phase transition, see 
Fig. [21 Within the antiferromagnetic phase, a first-order 
phase transition further occurs at smaller J/D, see the 
jump of the order parameter in Fig. [21 accompanied by 
a rearrangement of the Fermi surface. This is shown in 
Fig. [3l where we draw the quasiparticle (emission) spec- 
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FIG. 2: (Color online) The magnetic order parameter as func- 
tion of J/D for Tie = 0.92 (left panel) and — 0.7 (right 
panel). Notice that for ric = 0.92 the order parameter grows 
continuously below a critical J/D 0.6 - second order phase 
transition - until at J/D ~ 0.36 it jumps abruptly - first or- 
der transition. For ric ~ 0.7 only a first order transition with 
a jump from zero to a finite value of the order parameter is 
observed. 



We have solved the variational problem numerically us- 
ing, for numerical convenience, a flat conduction-electron 
density-of-states with half-bandwidth D, our unit of en- 
ergy. We do not expect that a more realistic density of 
states could qualitatively change the phase diagram that 
we find. Some technical details of the calculations are 
presented in the appendix, while here we just discuss the 
results. 



tral function at the chemical potential, defined by 

A{k)=- ldeA{k,e)^, (11) 

where /(e) is the Fermi distribution function at low tem- 
perature and 



A(k,6) = ^^^^e'''(^-^')5(i?„-i?o-e)(*o|(^cc.(R)cL + ^c/.(R)/L)l*«) 

RR' an ^ 

(vf,,! f z;,,(R') c^,^ + z:/,(R') f^A \^o) (12) 



r 



with V the number of lattice sites. |^o) and \'^n) are the 
ground state and an excited state, with energy Eq and 
En, respectively, of the variational Hamiltonian Ti.^, see 
Eq. (|A.21[) in the appendix. A(k, e) is calculated with 



a nearest-neighbor hopping on a two dimensional square 
lattice, though with variational parameters optimized us- 
ing a flat density of states at the same values of ric and 
J/D. 
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The k-points where ^(k, e) is large identify the effec- 
tive Fermi surface. We note that, in the paramagnetic 
phase, the Fermi surface is hole-Uke just as if the / spins 
do partecipate the Luttinger sum rules - two bands with 
1 + Tic < 2 electrons per site; one band empty and the 
other occupied by 1 < 1 + ric < 2 electrons. The same 
feature is also found beyond the second order phase tran- 
sition. However, for J/ D below the first order phase 
transition, the Fermi surface changes topology and be- 



come electron-like, as if the /-electrons disappear from 
the Fermi surface. Comparing the phase diagram Fig. [T] 
with the one obtained by VMO^i, we find that the two 
agree well, even quantitatively!^ In order to identify the 
origin of the Fermi surface rearrangement, it is conve- 
nient to write the general expression of the variational 
Hamiltonian 7i,, see Eq. (jA.2ip . of which j'l'o) is the 
ground state. In momentum space and within the mag- 
netic Brillouin zone 
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(13) 



where ek is the energy dispersion of the conduction elec- 
trons, 



"^k(T' /kcr' "^k+Qcr' /k+Q(T ) ' 



a Fermi spinor, its hermitean conjugate being "iAko-' Q 
Neel magnetic vector, and all the Hamiltonian parame- 
ters are variational but ek- 
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FIG. 3: (Color online) The conduction electron spectral func- 
tion at the chemical potential for a two-dimensional square 
lattice. Panels (a)-(b)-(c) show the evolution of the spectral 
function ^(k) at the chemical potential for ric = 0.92 in the 
paramagnetic phase, panel (a) with J/D = 0.8, right after the 
second-order transition, panel (b) with J/D = 0.4, and finally 
below the first-order transition, panel (c) with J jD = 0.16. 
Panels (d)-(e) show the same evolution with ric = 0.7 where 
there is only the first-order transition. 

In Fig. [5] we plot the variational bands in the antifer- 
romagnetic phase below and above the first order phase 
transition. In agreement with the interpretation given 
by Watanabe and Ogata in Ref. [2l|, the bands in the 
antiferromagnetic phase at low J jD can be sought as 
antiferromagnetically split c and / bands very weakly 
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FIG. 4: (Color online) Evolution of the band structure of the 
optimized variational Hamiltonian Eq. (I13|l for ric = 0.92 as 
a function of JjD and across the first order transition. From 
top left to bottom right panel: JjD = 0.1, JjD = 0.2, JjD = 
0.3 (below the first-order transition) and J/ D = 0.36 (above 
the first-order transition). 



hybridized, panels (a) and (b) in Fig. O while those at 
larger J / D as strongly hybridized c and / bands weakly 
antiferromagnetically split, panels (c) and (d) in Fig. O 
The main control parameter of the transition is the rel- 
ative strength of the /-orbital energy, e/ in (fTH]) . with 
respect to the antiferromagnetic splittings, mostly aAI 
in dni). 

Above a critical doping away from the compensated 
regime, we only find a single first-order phase transition 
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FIG. 5: (Color online) One-dimensional representation of 
the diflterent variational band structures in the two magnetic 
phases close to nc = 1, drawn in the magnetic Brillouin zone. 
Small J jD phase: panel (a) represents non-hybridized c and 
/ bands split by antiferromagnetism; panel (b) what happens 
once a small hybridization is switched on. Large J/_D phase: 
panel (c) represents non-magnetic hybridized c and / bands 
in the folded Brillouin zone; panel (d) what happens once a 
small antiferromagnetic order parameter is switched on. 

transition, see Fig. [21 directly from a paramagnet at large 
J/I?, with a band structure similar to panel (c) in Fig. O 
unfolded in the whole Brillouin zone, to an antifcrromag- 
net with a band structure similiar to panel (b) in Fig. [5] 
In other words, this phase transition is accompanied by 
a drastic reconstruction of the Fermi surface. 



IV. FERMI-SURFACE RECONSTRUCTION VS. 
MAGNETISM 

The variational phase diagram. Fig. [1] shows that the 
onset of magnetism is not necessarily accompanied by a 
Fermi surface reconstruction. Viceversa, one could spec- 
ulate that the latter might not require magnetism, which 
would be the case if the Fermi-surface change were caused 
by the /-electron localization^. This aspect makes worth 
investigating the properties of the variational wavefunc- 
tion ([2]) preventing antiferromagnetism, which amounts 
to assume Arn(R-) in Eq. ([3]) independent of R and j^I'o) 
a paramagnetic Slater determinant. 

At first sight, one would not expect to find anything 
special varying J jD in the paramagnetic sector. In fact, 
we previously mentioned that the change of the Fermi 
surface within the magnetic phase reflects essentially the 
change of the band structure, which, in turn, depends 
variationally only on the value of the /-orbital energy 
with respect to the magnetic splitting, respectively e/ 
and 2Af in Eq. (jl3p . Therefore, without magnetism, i.e. 
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FIG. 6: (Color online) Variational energy as function of J jD 
at ric = 0.8 in the paramagnetic sector. A kink is visible at 
J jD ~ 2.1. We note the finite curvature of the energy at low 
J/D, which, as we checked, is compatible with second order 
perturbation theory. 



(a) (b) 




FIG. 7: (Color online) Spectral function at the chemical po- 
tential for ric. = 0.8 and JjD above, panel (a), and below, 
panel (b), the critical value. 



M = 0, the topology of the band structure must remain 
invariant whatever JjD ^ is, as we indeed find. Never- 
theless, even in this case, we do observe a very weak first 
order phase transition for values of JjD slightly smaller 
than those at which the first order transition occurs when 
we allow for magnetism, as shown by the behavior of the 
variational energy in Fig. [S] Remarkably, at this transi- 
tion the momentum dependent spectral function of the 
conduction electrons at the chemical potential changes 
abruptly, see Fig. [T] For JjD above the critical value, 
the Fermi surface includes the /-electrons, while, below, 
it does not, exactly as we find when magnetism is present. 
However, this change occurs now not because the band 
structure is modified but because the spectral weight of 
the conduction electrons at the Fermi energy changes dis- 
continously. Indeed, looking carefully at the momentum 
distribution in Fig. [7^, one can distinguish two sheets of 
the Fermi surface, a small one, which corresponds to the 
non-interacting conduction electron Fermi surface, and a 
large one that includes also the / electrons. Across the 
transition, it is the relative weight of these two sheets 
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FIG. 8: (Color online) Behavior of the average Kondo ex- 
change, K/D, and hopping, T/D (in units of D, half the 
conduction bandwidth). 




FIG. 9: (Color online) Evolution of the uniform magnetiza- 
tion as function of an external magnetic field applied in the 
paramagnetic phase {J / D = 0.45 at ric = 0.88). Insets show 
the spectral functions for the majority (top panel) and minor- 
ity (bottom panel) spins across the metamagnetic transition. 



that change discontinuously. We believe that this must 
be regarded as a manifestation of an /-localization, or, 
better, of an orbital-selective localization, as proposed in 
Refs. [l3 andd^. This result also demonstrates that the 
rearrangement of the Fermi surface observed along the 
first-order line in the phase diagram Fig. [1] is caused by 
the /-electron orbital-selective localization rather than 
by magnetism. 

Inspection of the behavior of the average Kondo ex- 
change and hopping, Fig. [51 shows that the "localized" 
phase has a better conduction-electron hopping energy, 
while the "delocalized" one a better Kondo exchange. 
This suggests that the abrupt change of the Fermi sur- 
face is primarily consequence of the competition between 
the conduction electron band-energy and the Kondo ex- 
change, and not of the commonly invoked compction be- 
tween Kondo and RKKY interactions. 

In light of these results, also the transition lines in the 
phase diagram. Fig. [TJ assume a different meaning. The 
first-order line that separates the paramagnet from the 
antiferromagnet is primarily due to the /-localization, 
magnetism being just its by-product. On the contrary, 
the second-order line close to the compensated regime 
is more likely to be interpreted as a Stoner's instability 
of the paramagnetic Fermi-liquid, driven by the nesting 
property of the Fermi surface at ric ~ 1. Across this 
second-order phase transition, the Fermi surface changes, 
smoothly, following the spin splitting of the bands. 



V. METAMAGNETISM 

Another indirect signal of the /-localization can be 
found by studying the behavior of the paramagnet in 
the presence of a uniform magnetic field. Indeed, if the 
/-orbitals are close to a Mott localization, they are also 



very prompt to order magnetically. Let alone, they would 
prefer some magnetic order along with the structure of 
the RKKY exchange, in our bipartite lattice model not 
far from half-filling the natural candidate being a Neel or- 
dering. However, in the presence of a magnetic field, they 
could equally prefer to order ferromagnetically. In other 
words, it is plausible to foresee that the /-localization 
could be driven by a weak magnetic field, the weaker the 
closer the orbital-selective Mott transition is, thus accom- 
panied by a sharp increase of magnetization, so-called 
metamagnetism, as well as by a discontinuous change of 
the Fermi surface. 

This expectation is confirmed by our variational cal- 
culation. In Fig. [9] we show the evolution of the uniform 
magnetization as function of the applied magnetic field 
in the paramagnetic phase at J /D = 0.45 and ric = 0.88. 
Indeed, as function of the magnetic field, we do find a first 
order phase transition that is accompanied by a abrupt 
increase of the magnetization as well as by a discontin- 
uous change of the conduction electron Fermi surface, 
specifically of the majority spin one. In fact, since the 
critical field is smaller than the Kondo exchange J, once 
the / electrons localize and their spins align with the 
external field, the effective Zeman field felt by the con- 
duction electrons is opposite to the applied one. Conse- 
quently, the Fermi surface of the majority spin becomes 
smaller than the minority spin one, contrary to the case 
for external fields below the metamagnetic transition, 
which is what we find, although hardly visible in Fig. [9l 



VI. CONCLUSIONS 

We have calculated within the Gutzwillcr approxima- 
tion the phase diagram of the Kondo lattice model as 
function of the conduction electron density and of the 
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Kondo exchange J. The novel feature of our approach 
with respect to earher ones is that the Gutzwiller projec- 
tor acts on all the electronic configurations of each / or- 
bital plus the conduction state to which it is hybridized. 
This allows to include additional local correlations be- 
tween / and conduction electrons, specifically those that 
favour singlet pairing among them. Summarizing our 
variational results, we have found that: 

• there exists an orbital selective Mott localization 
of the / electrons accompanied by a discontinuous 
change of the Fermi surface; 

• away from any nesting instability, this first-order 
transition in accompanied by magnetism; 

• on the contrary, when the conduction electron 
Fermi surface is perfectly or almost perfectly 
nested, magnetism occurs before the /-localization, 
via a second order transition with a continuous 
change of the Fermi surface; 

• the /-electron Mott localization can be also induced 
by a uniform magnetic field, in which case it is re- 
vealed by a metamagnetic transition at which the 
magnetization jumps and the Fermi surface changes 
discontinuously. 

These findings bridge between the cluster dynamical 
mean field theory results of Refs . [iTtllsl and the vari- 
ational Monte Carlo ones of Ref. l2ll . and suggest that 
generically, i.e. without nesting, magnetism is a by- 
product of the /-electron Mott localization rather than 
the outcome of the competition between Kondo screen- 
ing and RKKY interaction. We must mention that the 
weak first-order character of the Mott transition that we 
find might be a spurious outcome of the variational pro- 
cedure, so that we can not exclude that in reality such a 
transition is continuous. 

The question we can not address, since ours is a 
variational approach for the ground state, concerns the 
anomalous thermodynamic behavior observed around the 
magnetic transition. In other words, we can not establish 
whether such a behavior is associated with the incipient 
magnetisroi^ or is just a consequence of the /-electron 

localization^iSiii or better of the orbital selective Mott 
localization.iSJLSOSdo 
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APPENDIX: THE GUTZWILLER 
APPROXIMATION IN DETAIL 

In this appendix we present some technical details of 
the method that we have employed, which simplify con- 
siderably all calculations. We start by the definition 



of the Gutzwiller wavefunction Eq. ([2]) with the general 
Gutzwiller operator Vr of Eq. ([3]). We assume that the 
average value of the local single-particle density-matrix 
operator, Cro- in Eq. ([7]), on the uncorrelated Slater de- 
terminant wavefunction ^'o is diagonal in terms of the 
operators d\^^ and ^2^^^, related by a unitary transfor- 
mation to the original ones, c^^^ and /r^.. In other words, 
for a,b — I, 2, 



(*0| dl-^^d^j^^ |*o) = Sab (5a/3 n^aRa, 



(A.l) 



where < n'^R^ < 1 are the eigenvalues of Cr^. We 
specify Vr to be of the form as in Eq. ^ , namely 



7'r = V Ar„(R)|r,R)(n,R|, 



Fn 



(A.2) 



where the states |r, R) are written in the original c-f 
basis and contain only one / electron while, by assump- 
tion, |n,R) are Fock states in the natural basis, namely 
in terms of di-d2- In other words, and dropping for sim- 
plicity the site-label R, a generic state \n) is identified by 
the occupation numbers Uaa = 0, 1, a = 1, 2 and a 
and has the explicit expression 

I-) - {^r {<r' i'lr {<r' 

We introduce the uncorrelated occupation-probability 
matrix P*^ with elements 



(^'o||m)(n||«'o) =<5„™P°, (A.3) 



where 



n n 



,0 



1-n' 



aa) 



(A.4) 



We also introduce the matrix representation of the oper- 
ators d^a- and d\^, namely 

daa {da„)nm = (^l^aalm) , 

dL (dL)„„ = {n\dL\m) = {{m\da,\n)y , 

and assume that the variational parameters Apn in 
Eq. (|A.2p are the elements of a matrix A. With the above 
definitions, the two conditions Eqs. ([5]) and ([6]) that we 
impose, and which allow for an analytical treatment in 
infinite-coordination lattices, become^^ 



(*o|-pt-p i^g) ^ Tr(pOAt A 



(A.5) 



(vE-ol Pt p dLrffc^ |*o) = Tr (P" At A dl,d,^ 



E ^nAipAp„(mMLd 



Trim 



"fol di^d,., \^o) ^ Sab Sc^pnl^ 



(A.6) 
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If Eqs. (jA.5|) and (jA.6[) are satisfied, then the average 
value of any local operator O in infinite-coordination lat- 
tices isM^ 

(^lOI*) = (*o|7'tC'-p|*o) 

= Tr (P" X^OX)^J2 Pn Air Orr' Ar,„, (A.7) 
Tirr' 

where O is a matrix with elements 

Orv = (r|0|r'). 

In the mixed original-natural basis representation, the 
proper definition of the Z-factors in Eqs. ([9]) and pO)) 
changes into 

(M'ol 4 V d,, 1*0) = Zci. (*o| did,, |*o) 
= ^ci. (A.8) 

+ Z,2. {^ald^d^J^o) 
= Z,2an°,. (A.9) 



In other words, when calculating the average hopping, 
the operator cj^ effectively transforms into 



4 ~^ Zcla dl + Zc2a C^L 



(A.IO) 



1. Explicit formulas and connection with 
slave-boson mean field theory 

To further simplify the calculation, we introduce a new 
matrix in the mixed basis representation 



with elements 



= AVPO, 



(A.ll) 



(A.12) 



As we shall see, 4>rn corresponds to the slave-boson 
saddle-point value within the multiband extension of 



the Kotliar-Ruckenstein mean-field scheme recently in- 
troduced by Lechermann and coworkers^, which they 
named rotationally invariant slave-boson formalism. By 
means of the definition (|A.11[) the first condition ([5]) be- 
comes 



Tr (0t = ^ 



At 



1, 



which coincides with the saddle-point value of Eq. (28) 
in Ref. i41|. The other condition, Eq. (jA.6[) . becomes 



Tr VPO, 



Trim 




P' 



6/3 I 



SabSaf3 



(A.13) 



which is apparently different from the saddle-point value 
of Eq. (29) in Ref. HJ, which reads 



Tr 



Fnm 



(A.14) 



Therefore the equivalence between the rotationally 
invariant slave-boson formalism and the multi-band 
Gutzwiller approximation is not so immediate as claimed 
recently by Biinemann and Gebhard4^ Indeed the proof 
given by these authors suffers by a flaw^, eventhough, 
as we are going to show, the final conclusion seems to be 
right. In fact, we note that the two Fock states |n) and 
|m) in Eq. (jA.13p differ only because |n) has the orbital 
b with spin /? occupied but orbital a with spin a empty, 
while it is viceversa for |m), so that 




I n'b, (1 ~ nL) 



hence Eq. (jA.13p is actually equal to 



(1 - ng^) n„„ ^^^^^ 



Because of SatSap on the r.h.s this equation is equivalent equivalent to impose 
to Eq. (|A.14p provided ng^ ^ and 1. Therefore, 

if the average value of the single-particle density matrix ^Pj. ^^t 0^ = 1 (A 15) 

on the uncorrelated wavefunction I^Pq) has eigenvalues 

neither nor 1, the conditions Eqs. (jA.Sp and (|A.6p are Tr (j) d^^d^^j = (*ol dl^d^p l*o) 
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= SabSapn'l^, (A. 16) which coincide with Eq. (47) in Ref. fill at the saddle 

point. Finally we need to evaluate Zda and Zc2ct of 
which indeed coincide with Eqs. (28) and (29) in Ref.|4l| Eqs. fO|l and (jAjJ. We find that, for a = 1, 2, 
specialized to the natural orbital basis and evaluated at 
the saddle point. 

In terms of (p, the average of the local operator O, 
Eq. (jA.7|) . becomes 

(^-lOI*) = Tr(<?!)tO0) , (A.17) 



■ Tr \/po -/-t 4 . 



-,0 '^aa 



^ E \/S<^lr(r|4|rVr'™(m|d,Jn) 



(1 - nO J p^ 



^ 0lp(r|4|r')0p,„(m|d,Jn), 



(A.18) 



rr'nm 



which are analogous to those proposed by Lechermann 
and coworkers^, as shown by Biinemann and Gebhard.— 

2. The variational energy 

The average value of the Hamiltonian ([T|) is the sum of 
two terms, the average of the hopping Ho plus that of the 
Kondo exchange TCj. The latter is a purely local term, 
hence, provided Eqs. (jA.lSP and (jA.16P are verfified and 
in infinite-coordination lattices , can be written through 
(|A.17[) as (we reintroduce the site label R) 

(vf-l 7^, I*) = J ^ Tr |^0t(R) s^j, . S,R ■ (A. 19) 
I 



under the constraint that the single particle density ma- 
trix has eigenvalues n'^nc,- One readily realizes that |^'o) 

I 



where the Lagrange multipliers /Liabo-(R) are those that 
maximize the ground state energy, so that 

(*|7^o|*> = (*o|W*|*o). (A.22) 



A way to proceed is to use as variational parameters the 
matrices 0(R) and the n^-fi^^s,, which are related to each 
other by the conditions Eqs. (|A.15|) and (jA.16p . Then, 
the Slater determinant 1^*0) must be the one that mini- 
mizes the average value of the hopping, which is, through 
Eqs. (|XTO)) and ^EM, 



(A.20) 



that fuUfiUs such a property is actually the ground state 
of the following variational Hamiltonian 



(A.21) 



Once this task has been accomplished, one needs to min- 
imize the variational energy per site 

Evar = 4(*0|W* |*0) 



I 



<RR'>cr a, 6=1 



<RR'>CT a, 6=1 \ / 



EEE 

R a, 6=1 o- 



Ma6a(R) ( ^RaC^bRa " '^"6 "aRa I + H.C 
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+ J ^ Tr (^(/.t (R) S/R • S,R 0(R) j ,( A.23) 

with respect to 0(R) and n^n^ that fulfill Eqs. (lA.lsp 
and (|XT6)) . 

3. Technical remarks 

In order to parametrize the variational matrix ^i(R) 
(in what follows we assume a generic multiband Hamil- 
tonian) one can introduce a local Hamiltonian hn that 
acts on all possible local electronic configuration, and de- 
fine 

</.(R)V(R) = -— , (A.24) 

where 

= Tr (e-^''-) , 

is the local partition function and 1//3 a fictitious tem- 
perature. With this definition, the condition Eq. IA.16I 
becomes 

Cabin) ^ ^ Tr (e-'^'^- cl^c,^) ^ (vfolc^c.^l^o), 

^ (A.25) 
where a and b label both spin and orbitals. Therefore, 
the zero-temperature average value of the single-particle 



density matrix on Ivfo) must coincide with the its thermal 
average with the local Hamiltonian /ir. Given /ir, one 
calculates Ca6(R), Eq. (jA.25p . which can be diagonalized 
providing the definition of the natural basis: 

^ Tr (e-^''- di^d,^) = 6ab Uan- (A.26) 

In terms of /ir 

m) = Un^^, (A.27) 
V "R 

with a unitary matrix. The expressions of the Z 
renormalization factors are then obtained through 

^ Tr (^e^^'- Ui ci^ e"^''- d,^ 

= ZabiR) V"bR(l-"bR)- (A.28) 

We found that it is more convenient to use as variational 
parameters instead of the matrix elements of 0(R) those 
of the local Hamiltonian /ir and of the unitary matrix 
L/r. In the case of a paramagnetic wavefunction that 
does not break translationally symmetry, /ir and C/r are 
independent of R. On the contrary, for antiferromagnetic 
wavefunctions on a bipartite lattice, going from one sub- 
lattice to the other the role of spin f (j) is interchanged 
with that of spin I ("f). 
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